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We provide an analytical solution to the problem of scattering of electromagnetic radiation by 
a square- wave grating with a flat graphene sheet on top. We show that for deep groves there is 
a strong plasmonic response with light absorption in the graphene sheet reaching more than 45%, 
due to the excitation of surface plasmon-polaritons. The case of grating with a graphene sheet 
presenting an induced periodic modulation of the conductivity is also discussed. 
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I. INTRODUCTION 

Plasmonic effects in graphene is currently an ac- 
tive research topic. The strong plasmonic response of 
graphene at room temperature is tied up to its optical 
response, which can be controlled externally in differ- 
ent ways. The unique features of the optical conduc- 
tivity of single-layer graphene stem from the Dirac-like 
nature of quasi-particlesP^ and have been ext ensively 
studied in the past few years, both theoreticall ^ 1 | 2 | 4 tt ^ 
and experimentallyP^HH] including in the terahertz (THz) 
spectral rangeP^H 

Indeed graphene holds many promises for cutting edge 
THz applications, 29 which would be able to fill the so 
called THz gap. More recently the interest has been fo- 
cused to how graphene interacts with electromagnetic ra- 
diation in the THzP^^ One of the goals is to enhance 
the absorption of graphene for the development of more 
efficient photodetectors in that spectral range. This can 
be done in several different ways, by (i) producing micro- 
disks of graphene on a layered structure} 3 -^ (ii) exploit- 
ing the physics of quantum dots and metallic arrays on 
graphene) 32 1 37 1 (iii) using a graphene based gra tingp SEH 

(iv) putting graphene inside an optical cavityHlESl anc [ 

(v) depositing graphene on a photonic crystal.^ In cases 
(i), (ii), and (iii) the excitation of plasmons^^ is re- 
sponsible for the enhancement of the absorption. In case 
(iv) photons undergo many round trips inside the cavity 
enhancing the chances of being absorbed by graphene. In 
case (v) the authors consider a photonic crystal made of 
Si02/Si. In the visible range of the spectrum the dielec- 
tric constants of Si02 and Si differ by more than one or- 
der of magnitude and choosing the width of the Si02/Si 
appropriately it is possible to induce a large photonic 
band gap in the visible range. Combining the presence 
of the band gap with an initial spacer layer the absorption 
can be enhanced by a factor of four. In the case studied in 
Ref. [15] the optical conductivity of graphene is controlled 
by interband transitions. Although this workP^ focused 
on the visible spectrum, there is a priori no reason why 
the same principle cannot be extended to the THz. 



The physics of surface plasmon-polaritons in graphene 
has also been explored for the development of devices 
for optoelectronic applications! 4 ^* 5 -^ Such devices include 
optical switches^ and polarizers.^ It has been shown 
that metallic single- wall carbon nanotubes, which have a 
linear spectrum close to zero energy, can act as polariz- 
ers as well. 27 Theoretical studies of the optical response 
of graphene under intense THz radiation has also been 
performed.^^ 

In general, the conductivity of graphene is a sum of 
two contributions: (i) a Drude-type term, describing in- 
traband processes and (ii) a term describing interband 
transitions. At zero temperature the optical co nductiv- 
ity has a simple analytical expression! 1 * 2 * 4 ! 5 * 8 -! 5 - 7 -] In what 
concerns our study, the physics of the system is domi- 
nated by the intraband contribution,^ which reads 
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where, cjq = 7re 2 /(2h), T is the relaxation rate, Ep > is 
the Fermi level position with respect to the Dirac point, 
and uj is the radiation frequency. It should be noted that 
(Jd has a strong frequency dependence and is responsible 
for the optical behaviour of graphene in the THz spectral 
range. 

The problem we consider in this work is the scatter- 
ing of electromagnetic radiation by graphene deposited 
on a grating with a square- wave profile, as illustrated 
in Fig. [I] We further assume that the radiation is 
p— polarized (TM wave), that is, we take B = (0,5^,0) 
and E = (E x ,0,E z ). The square- wave profile is rather 
pathologic due to the infinite derivative at the edges of 
the steps and therefore the usual methods fail! 38 * 40 * 59 * 60 * 
An alternative route is to obtain the exact eigenfunc- 
tions in the region of the grooves.^ Fortunately, the ge- 
ometry of the problem is equivalent to that of a Kronig- 
Penney model appearing in the band theory of solids and 
for which an exact solution exists. We will show that 
Maxwell's equations can be put in a form equivalent to 
the Schrodinger equation for the Kronig-Penney model, 
and hence an exact solution for the fields in the region of 
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FIG. 1: Geometry of the problem. The incident angle is 
#2 and the incoming wave number is ko. The radiation 
is p— polarized. A top gate for the electrostatic doping of 
graphene can be arranged as a transparent electrode placed 
at some distance above the graphene sheet. Alternatively, a 
bottom gate can be placed below the dielectric substrate. 



the grooves is possible. 



In this work we assume that graphene is doped. In a 
practical implementa tion, t his can be achieved via gating 
or by chemical meansP^ES Remark that in a bottom gate 
structure the conductivity of graphene becomes position 
dependent along the x— direction!^ In this case, the anal- 



ysis presented below is incomplete (see, however, Sec. VI 
where a simplified model is analysed). Nevertheless, the 
results of Ref. 38 show that a position dependent conduc- 
tivity alone already induces surface plasmon-polaritons 
on graphene. Therefore, a spatially modulated optical re- 
sponse, when combined with a grating, will enhance the 
excitation of surface plasmon-polaritons, as confirmed by 
the analysis given in Sec. VI For the top gate configu- 



ration and for chemical doping, the problem of a spatial 
dependent conductivity does not arise!^ 



It can be argued that in the geometry we are consider- 
ing (see Fig. [I]) the portion of graphene over the grooves 
will be strained. Clearly, experimental setups can be 
designed as to overcome (at least partially) such effects 
e.g., by filling the grooves with a different dielectric and 
using a top gate or a chemically doped graphene. How- 
ever, in the bottom gate configuration, even with filled 
grooves, graphene will have a position dependent con- 
ductivity. A complete description should thus include all 
contributions: the grating effect itself, inhomogeneous 
doping, and the strain fields in graphene. In this work 
we assume that the graphene sheet remains flat. First, 
the optical conductivity will be taken homogeneous in 
order to study solely the effect of dielectric grating in 
the presence of uniform graphene. The additional ef- 
fect of periodic modulation of the graphene conductivity, 
with the same period as the grating, will be considered in 
Sec. EH We believe that our calculations are accurate for 
the cases of either electrostatic doping by a top gatep^ 
or doping by chemical means P^EsD 



II. EXACT EIGENMODES IN THE GRATING 
REGION 

For TM polarized light with angular frequency uo (see 
Fig. [T]) the Helmoltz equation assumes the simple form 

(Ay + no6uo 2 )B y (x,z) = 0, (2) 

where /jlq is the vacuum permeability and An = d 2 + d 2 . 
Since the boundary between the vacuum and the dielec- 
tric is piecewise, the Helmholtz equation holds true in 
the regions of width rd and (l—r)d with the appropriate 
change of the dielectric constant, e in the former and eo 
in the latter. For z > h the dielectric constant is homo- 
geneous and equal to eo, whereas for z < it is e = eoe r , 
and e r is the relative permittivity. 

We search solutions of Eq. [2] in the form 



Bo, 



B y (x, z) = X(x)e 



±iAz 



(3) 



where A is a constant. With this ansatz, the Helmholtz 
equation reduces to 

d 2 x X(x) = (A - fi euj 2 )X(x) , 

whose solution is 



(4) 



X(x) 



and where k" 2 - is given by 



fc 2 = /xoejw 2 - A 2 , 



(5) 



(6) 



with the subscript j = 1 [j = 2] referring to the region 
(1 — r)d [rd]. Putting all pieces together, the field com- 
ponent By in the regions of the grooves reads 



= [Aje 



■Bje 



(7) 

The determination of A in Eq. ([7]) leads to an eigenvalue 
problem that will be considered in the next section. We 
note that A can be either real or pure imaginary: re- 
als values correspond to propagating diffracted orders, 
whereas imaginary ones correspond to evanescent waves. 
From Maxwell's equations it follows that the electric field 
components in the region < z < d read: 



[A je 



Bje 
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-ikjX~\^±iAz 



-ikjX~\^±iAz 



(8) 



(9) 



As shown in what follows, there is a relation between 
the coefficients Aj and B^ which is obtained from the 
solution of the eigenvalue problem. 



III. TRANSFER MATRIX AND THE 
EIGENVALUE EQUATION 

Along the x— direction, and for < z < h, we have a 
stratified medium with alternating dielectric constants. 
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Then we can relate the amplitudes Aj and Bj with A/+i 
and Bj+i using the boundary conditions at the interfaces 
of the two dielectrics. Furthermore, using the Bloch theo- 
rem we can find an eigenvalue equation for the parameter 
A. The amplitudes in the different regions of the strat- 
ified medium are represented in Fig. [2] The continuity 



ei = eo £2 = e e r 
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FIG. 2: Effective stratified medium along the x— direction in 
the region < z < h. 



Combining Eqs. (12) and (18), we arrive at the following 
result 



A 2 
B 2 



^(k 2 rd)K- 1 K 1 ^[k 1 d(l - r)]K~ 1 K 2 



A 4 
B 4 



(19) 



where a global phase was absorbed in the coefficients Ai 
and B4. The transfer matrix, 



T = [QihrdjK^K^lhdil - r)\K^K 2 ] 



(20) 



propagates the field amplitudes through the ID crystal 
in the x— direction, 



of the ^-component of the electric field at the boundary 
x = rd imposes a relation between the amplitudes A and 
B in regions j = 2, 3 (refer to Fig. [2] for the definition of 
the different regions), i.e., 

— [A 2 e ik2rd - B 2 e~ ik2rd ] = /ci[A 3 e i/cir(i - ^^~ ik ^ 



B 3 e~ 



(10) 

Note that according to the definition in Eq. (|6| , we must 
have ks = k\ and k^ = k 2 . At the same boundary, the 
continuity of the magnetic field yields 



A 2 e ik2rd + B 2 e~ ik2rd = A 3 e ikird + B 3 e~ ikird . 



(11) 



It is convenient to write these two equations in matrix 
form as follows 



A 2 
B 2 



where 



= ^(k 2 rd)K^ 1 K 1 ^(-k 1 rd) 



K 2 



A3 
B 3 
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(12) 

(13) 
(14) 
(15) 



Similar continuity conditions apply to the boundary at 
x — resulting in the following constraints 



k x [A 3 e ikld - B 3 e- ikld ] = — [A A e ik2d - B A e~ ik2d ) , (16) 

6 r 



A 3 e ikld + B 3 e- lKia = A±e 



ik\d 



B 4 e 



(17) 



As before, these equations can be written in matrix form 

as 



Az 
B 3 



^(k^K^K^i-hd) 



Ai 



(18) 



Ai 
Bi 



A 2 
B 2 



(21) 



[Note that det (T) = 1.] On the other hand, by virtue of 
the Bloch theorem we have: 



Ai 
Bi 



O ikodsin i 



A 2 
B 2 



(22) 



where ko sin 0i is the x— component of the wave vector of 
the incoming electromagnetic radiation (Fig.JT]). We thus 
arrive at the important intermediate result 



O ikod sin t 



A 2 
B 2 



A 2 
B 2 



(23) 



The compatibility condition of Eq. (23) provides the 
eigenvalue equation 



2cos(/c dsin<9 i ) = Tr(T) . 



(24) 



or explicitly, 



2 cos (/cousin 0, 
f 2 k 2 



2cos[kid(l — r)] cos(k 2 dr) 



k\ 



e r kik 2 



sin[kid(l — r)] sin(k 2 dr) . 



(25) 



Eq. (25) allows for the determination of the permitted 
values of A 2 and is very similar to the eigenvalue equation 
of the Kronig-Penney model of electron bands. 

The relation (12) allows us to express the coefficients 



A 3 and B3 in terms of A 2 and B 2 . Thus, the function 
X(x) over the whole unit cell can be written in terms of 
A 2 and B 2 only. Furthermore, Eq. (23) gives a relation 
between the coefficient A 2 and B 2 , 



B 2 



O ikodsin 9i 



'12 



(26) 



where T\ 2 and T 22 are two matrix elements of the transfer 
matrix T. Therefore, the function X(x) over the unit 
cell is proportional to the only coefficient, A 2 . Note that 
the matrix elements contain A 2 , therefore we shall label 
the possible functions X(x) by index £ running over all 
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possible eigenvalues (±A^). The function Xg{x) in a unit 
cell has the form 



*) = { 



A 2 e ik2X + B 2 e- ik2X , < x < rd , 
A 3 e iklX + B 3 e~ iklX , rd<x<d. 



(27) 



Equation (12) can be writtent explicitly as 



where 



A 3 = aA 2 + bB 2 , 
B 3 = cA 2 +fB 2 , 



- i{ - kl - k ^ rd {k 2 + k 1 e r ), 



f 
b 
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(-fe 2 + feic r ) . 



i(k\+k2)rd 



( — fc 2 + fcie r ) . 



Since the wave numbers fcj can be complex, / is not nec- 
essarily the complex conjugate of a; the same applies to 
b and c. In terms of these coefficients, X £ (x) reads: 



A 2 e ik2X + B 2 e~ ik2X , < x < rd , 

(aA 2 + 6B 2 )e ifcias + (28) 
(cA 2 + fB 2 )e~ iklX , rd<x<d. 



X £ (x) 



We also note that Eq. ([26]) allows to replace £? 2 in Eq. 
(28). Finally, the magnetic field in the region < z < h 



(hereafter denoted as region II) has the form, 

= ^X £ (x)(C l e iA ' z + D £ e- iA ' z ), (29) 
i 

where the summation is over all the eigenvalues deter- 
mined from the solution of Eq. (25) and Ci and Dg are 



some coefficients that will be determined in the next sec- 
tion. With this we conclude the exact solution for the 
eigenmodes in the grating region. 



IV. SOLUTION OF THE SCATTERING 
PROBLEM 

We now derive the equations for the scattering problem 
represented in Fig. [I] For z > h (region I) the magnetic 
field is written as 



n— — oo 

and for z < (region III) we have 



, (30) 



x-p<?>z) 



(31) 



where 



= ko sin 0i — 27rn/d, 



v2 • 



< a n 



(32) 
(33) 



with q p — uj/vp and the definition qi = ko; n is an inte- 
ger, n G [—00,00]. Since we have four sets of unknown 
amplitudes, Ce, Di, R ni and T n , and four boundary con- 
ditions, two at z = h and other two at z = 0, we have 
a linear system of equations that can be solved in closed 
form if truncated to some finite order, iVg, which is the 
number of the eigenvalues needed in Eq. ( 29 ) for an ac- 
curate description of By 1 (we typically used Nz ~ 20). 
The boundary conditions at z = are 



B"\x,z 



0) 
0) 



E I J(x,z = 0), 



B I y I (x,z = 0), 



whereas those at z = h read 



E I x (x,z = h) = E I J(x,z = h), 
B I y {x,z = h) - B I y I {x,z = h) = 

- ^oa D E"(x,z = h) . 



(34) 
(35) 



(36) 



(37) 



The latter represents the magnetic field discontinuity 
across the graphene sheet 

These boundary conditions are x dependent. Since the 
system has period d, we can eliminate this dependence 
by multiplying the boundary conditions by e~ lOLrnX and 
integrating over the unit cell. After some algebra, we 
arrive at 



Y.^C. + f^Dd = 



E' 



(-); 



(38) 
(39) 



where 



r(±) 



a {±) 



with 
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X£r, 



1 f d 

= - J dxX £ (x)e 



dx^e-^*. 
e(x) 



(40) 
(41) 



Here e(x) is defined as e [eo] for x = rd [d(l— r)]. Eqs. (38) 
and ( 39 ) form a linear system of equations for the ampli- 



n=— 00 



tudes Ci and from which the reflectance and trans- 
mittance can be computed. Taking Ni odd, the integer m 



5 



belongs to the interval m e [-(N £ -l)/2, (N £ -l)/2\. The 
transmittance and the reflectance amplitudes are given 

by 



d=10nm, e=4 



=0.25, 0=0 E F =0.2 eV, r=1 meV 



(42) 



x ^2eoSl tm -fi(Cte iA ' h - D £ e- iA * h ) . (43) 

Since Xg{x) is a sum of exponentials, the functions xim 
and can be determined in closed form, which saves 
computational power. Explicit equations for xtm and 
Qtm are given in the Appendix [Al 



V. RESULTS FOR HOMOGENEOUS 
GRAPHENE 

We now provide a number of results obtained from the 
solution of Eqs. (38) and (39). In Fig. [3] we depict the 
absorbance (a), reflectance (b), and transmittance (c) at 
normal incidence. For the parameters considered only 
the specular order exists; all the other diffraction orders 
are evanescent. Then, the absorbance is defined as 



A=1-\R \ 



ft 



(2) 



(1) 



|r | 



(44) 



A resonance is clearly seen in the absorbance curves at 
a frequency of about 16 meV, corresponding to the ex- 
citation of a surface plasmon-polariton. The position of 
the resonance depends on several parameters, given in 
the caption of Fig. [3] and chosen as typical values ap- 
propriate for THz physics. The thin dashed curve cor- 
responds to graphene on a homogeneous dielectric (no 
grating, r = 1). The effect of increasing the depth of the 
grooves, from h = d/10 up to h = d is to produce an en- 
hancement of the absorption, which for h = d is almost 
of 25% (in this case we are considering r = 0.25). 

The dispersion of a surface plasmon-polaritons in 
graphene, when the sheet is sandwiched between two 
semi-infinite dielectrics, is, in the electrostatic limit, 
given b}E^ 



2a^ Epchq , 



(45) 



where 



Airehc 



and e = eo(l + e r )/2. Taking q = 2ir/d, the smallest 
lattice wavevector, Eq. (45) predicts, for the parameters 



used in our calculation, a plasmon-polariton energy of 



hoj = 12.1 meV, 



(47) 
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FIG. 3: Absorbance (a), reflectance (b), and transmittance 
(c) at normal incidence. The parameters are Ef — 0.2 eV, 
d=10 /mi, r = 0.25, 0i = 0, e r = 4 (corresponding to SiCb), 
and r = 1 meV. The thin dashed line is the case r — 1. For 
reference, we note that 1 THz corresponds to an energy of 
4.1 meV. 



which is in the ballpark of the numeric result. We should 
however note that the position of the resonance depends 
also on the parameter r, which is not captured by the 
above formula for the dispersion. We can include the 
effect of the parameter r using an interpolative formula. 
We define a new e as 



(2 - r + e r r) . 



Using this formula we obtain 



Huj = 16.3 meV , 



(48) 



(49) 



a much better approximation to the numerical value (15.9 
meV). To obtain the exact frequency of the resonance 
we have to compute the plasmonic band structure due 
to the periodic dielectric, that is, the band structure of 
a polaritonic crystal.^ In the Appendix [b] we give the 
polaritonic band structure of graphene on a square-wave 
grating. 

In Fig. [4] we depict the absorbance (a), reflectance (b), 
and transmittance (c), as function of frequency, for dif- 
ferent values of r, from r = 0.1uptor = 0.6, and keeping 
d = h, that is the limit of deep grooves. As r increases the 
position of the resonance shifts to the left. This happens 
because the width of the dielectric underneath graphene 
is increasing with r. Then, according to Eq. (46), the ef- 



(46) fective dielectric constant of the system increases, follow- 



energies. 



ing from Eq. (45) that the resonance shifts toward lower 

in t: 



le case of 



0.6 two resonances are seen 



in the frequency window considered. They correspond 
to the excitation of surface plasmon-polaritons of wave 
numbers 2ir/d and An/d. Then, according to Eq. ( |45| ), 
the position of the second resonance should be y/2 times 
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FIG. 4: Absorbance (a), reflectance (b), and transmittance 
(c) at normal incidence. The parameters are Ef — 0.2 eV, 
d=10 fjbm, h = d, 0i = 0, e r = 4, and r = 1 meV. The thin 
dashed line is the case r = 1. 



the first resonance frequency. From the figure, the en- 
ergy of the first resonance is fvjj = 13.6 meV, while that 
of the second one is Huj = 19.3 meV; now we note that 
y/2 x 13.6 = 19.2 meV, in agreement with the numerical 
result. We also note that the absorption of graphene is 
largest when r = 0.5, the symmetric case (1 — d = 0.5), 
reaching a value higher than 30%. Although we do not 
show it here, we have found that the absorption grows 
monotonically with increasing h and attains a maximum 
at a value larger than 45% for h = 2d. 

Finally, we note that the reflectance curves have Fano- 
type shaped This is due to the coupling of the exter- 
nal radiation field with the excitation of Bragg modes 
of surface plasmon-polaritons in periodically modulated 
structures. 

An important quantity when dealing with plasmonic 
effects is the enhancement of the electromagnetic field 
close to the interface of the metal (in this case graphene) 
and the dielectric (the square- wave grating). The spa- 
tial intensity of the diffracted electromagnetic field from 
graphene is depicted in Figj5j For the case off-resonance 
[Figj5ja)] the magnetic field has a relatively small am- 
plitude and is almost homogeneous along x-axis, while 
along z-axis exhibits a discontinuity at the plane z = 
occupied by the graphene layer, in accordance with the 
boundary conditions. 

The situation changes dramatically when the incident 
wave frequency coincides with that of a surface plasmon- 
polariton resonance [Figsjgb-d)]. In this case the am- 
plitude of the electromagnetic field in the vicinity of 
the graphene layer drastically increases. The electro- 
magnetic field amplitude is maximal when the surface 
plasmon-polariton resonance occurs for the first har- 
monic [Figj5jb)], that is for q = 27r/d, while it gradu- 
ally decreases as the harmonic's order increases (compare 
with Fig|5jc) for second and Figj5|d) for third harmon- 



FIG. 5: Spatial dependence of the magnetic field \B y (x, z)\ 2 , 
diffracted from the graphene layer on top of the square-wave 
grating. We have considered a groove depth h = d/2 and 
different frequencies of the incident wave: huo = 13meV (a), 
huj 16meV (b), Huj 23meV (c) and hw 28.4 meV (d). 
The other parameters of the structure are the same as in Fig. 

El 



ics, respectively). It is this enhancement of the electro- 
magnetic field due to the excitation of surface plamons- 
polaritons in the vicinity of a metallic interface that is at 
the heart of sensing devices. 



VI. MODULATED DOPING: A SIMPLE MODEL 

As discussed in the Introduction, doping graphene us- 
ing a bottom gate when the material lies on a grating 
leads to a position dependent conductivity, which is pe- 
riodic along the x— direction. In this section we consider 
a simple model where, in addition to the grating itself, 
the conductivity is position dependent and periodic. In 
our model the conductivity is given b}E^ 

2ttx 

<j(x) = <td(1 — k sin — ^— ) , (50) 



d 



where k is a parameter controlling the degree of inhomo- 
geneity. We also take r = 0.5 in our figures. This partic- 
ular choice of r and a{x) makes the latter commensurate 
with the profile of the grating. From the point of view 
of the calculation we have to replace gb by a(x) in the 
formulas given above. The conductivity is then expanded 
in Fourier series as 



a i27vxp/ d 



where 



1 f d 

- / dx<j(x)t 
d Jo 



-i27rpx/d 



(51) 



(52) 



For the case of the profile defined in Eq. ( 50 ) the Fourier 



series reduces to three terms only, those referring to 
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p = 0, ±1. A non- sinusoidal profile for a(x) will have 
more harmonics than just these three. Therefore, our 
model can also be considered the first term in the Fourier 
expansion of a more complex profile for <j(x). 

Relatively to the previous case of homogeneous con- 
ductivity, what changes in the equations is the form of 



the function 



9i 



(±) 



which now reads 



dim I X£m T — ^mAi 



Teo^^) 6 ™' ( 53 ) 



where 



— ^£ra — + — ^£ m _i . (54) 



The results for the absorbance, reflectance, and trans- 



d=10nm, E F =0.2 eV 



r=1 meV,e=4,e=0 
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FIG. 6: Absorbance (a), reflectance (b), and transmittance 
(c) at normal incidence. Three distinct situations are consid- 
ered: inhomogeneous conductivity, with grating (solid line) 
and without grating (dashed line), and homogeneous conduc- 
tivity with grating (dotted-dashed line). Other parameters as 
in Fig. [4] 

mittance are given in Fig. [6| considering the case where 
k = 0.5, a fairly large value, corresponding to deep 
groves. The dashed line is the case where graphene has 
an inhomogeneous conductivity on top of a homogeneous 
dielectric (i.e., no grating). (This situation is artificial 
and is only included for the sake of comparison.) Two 
resonant peaks can be seen, with the second peak be- 
ing much smaller than the one at lower energies. The 
dashed-dotted line is the case where we have homoge- 
neous graphene on the grating, and the solid line cor- 
responds to the case where we have graphene with an 
inhomogeneous conductivity on the grating. The main 
effect is a shift of the position of the resonant peaks to- 
ward lower energies and an enhancement of the absorp- 
tion, which bears its origin on the combined effect of the 
grating and of the periodic modulation of the conductiv- 
ity. In the absorbance panel the energy of the resonance 



labelled (2) is y/2 larger than that of the resonance (1). 
Thus, the first resonance corresponds to an excitation of 
a surface plasmon-polariton of wave number q = 2ir/d 
whereas the second peak corresponds to the excitation 
of surface plasmon-polariton of wave number q = 4rr/d. 
When we compare the dashed curve to the solid one we 
also note the prominence of the second plasmonic peak 
in the latter case, which is almost as intense as the low 
energy one. This is a consequence of the square profile 
and the dielectric nature of the grating. In metallic grat- 
ings the permittivity has a strong frequency dependence 
favouring the the occurence of the most intense plasmonic 
resonances at lower energies 



VII. CONCLUSIONS AND FUTURE WORK 

We have shown that a flat sheet of graphene on top of 
a square-wave grating exhibits strong plasmonic behav- 
ior. With the help of the grating one can create a sur- 
face plasmon-polariton resonance. The effect is more pro- 
nounced in the case of deep grooves, where absorbances 
higher than 45% are attainable. We note that the pa- 
rameter d is easy to control experimentally, providing a 
convenient way of tuning the position of the plasmonic 
resonance. The reflectance curves show a Fano-type line 
shape, which is manifest of the coupling of the exter- 
nal electromagnetic field (with a continuum of modes) to 
the surface plasmons in graphene (occupying a relatively 
narrow spectral band). The calculations we have per- 
formed assume that graphene is homogeneously doped 
(excluding Sec. VI), which is a crude approximation for 
the case of a bottom gate configuration of graphene dop- 
ing. Therefore, our results are only directly applicable to 
the cases of either chemical doping or doping by a top 
gate. This type of gating is within the state-of-the-art 
In the particular case of bottom gate, we have to con- 
sider both the effect of inhomogeneous doping and the 
effect of strain (if the grooves are not filled with another 
dielectric). Both effects result in a position dependent 
conductivity. Then, the calculation of the properties of 
surface plasmon-polaritons requires the evaluation of the 
doping profile and the strain field. Once these are deter- 
mined, one can use a Fourier expansion of the graphene 
conductivity as we have considered in our phenomeno- 
logical model. As these preliminary results show, the 
coupling of the external wave to the surface plasmon- 
polaritons can be significantly enhanced in the presence 
of both the dielectric grating and periodic modulation 
of the conductivity. Detailed calculations for realistic 
graphene conductivity profiles are the goal of a future 
work. 
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Appendix A: Explicit form of Xin 



The calculation of the integral in Eq. (40) has to be 



divided into two pieces. The function xtm can be written 

as 



_ (!) , ( 2 ) 
Xim — Xim "r Xim i 



(Al) 



where 



XI 



(i) 



d{k 1 - a m ) 
A 2 c 



d(ki + a m ) 
. B 2 b 
d(h - a m ) 



B 2 f 



d(ki + a m ) 



-i(k 1 +a m )d _ -i(k 1 +a rr i)rd 



and 



(2) 



+ i 



d(k 2 - a m ) 
B 2 



i(k 2 -ctm)dr 



i) 



d(k 2 + a m ) 



-i(k 2 +a m )d' 



r - 1) . (A3) 



Since the square profile has a piecewise structure, the 
function f|g m is obtained directly from the function xim 

as 



^irn — 



(1) (2) 
m 

eo e r eo 



(A4) 



Appendix B: Polaritonic spectrum 

The polaritonic spectrum, hu(q), of the surface 
plasmons-polaritons (SPP) in graphene on a square-wave 
grating is represented in Figj7| where q = ko sin 0^ is the 
Bloch wavenumber. Here, and in order to avoid the ap- 
pearance of an imaginary part of the eigenvalues, corre- 
sponding to surface plasmon-polariton damping, we con- 
sidered graphene without disorder (r = 0). As expected, 
the periodicity of the grating induces a band structure 
in the SPP spectrum of graphene, showing energy gaps. 
At the same time, for normal incidence (q = 0) the fre- 
quency of the second band (Hlj w 15.85 meV) almost co- 
incides with the numerically obtained resonant frequency 
hcu ~ 15.9 meV. A similar good agreement happens for 
the frequency of fifth band and that of second resonance 
(hoj w 23.1 meV and Huj ~ 23meV, respectively), as well 
as for frequency of sixth band and that of third resonance 
(hu « 28.3 meV and hu w 28.4 meV, respectively). 




q, jLim 



FIG. 7: Spectrum of surface plasmons-polaritons in graphene 
on a square- wave grating. We have considered Ef = 0.2 eV 
and r = 0. The other parameters are: d = 10 /xm, r = 0.25, 
e r — 4. The gray region represents the "light cone", limited 
by the light lines uj — c\q\ (dash-doted lines), inside which the 
excitation of surface plasmons-polaritons by incident light is 
possible. The frequencies, corresponding to the numerically 
obtained resonant minima of the transmittance for normal 
incidence, 0i = 0, are represented by black circles on the line 
q = 0. 
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